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Abstract 

The purpose of this article is to give the intelligible procedure for undergraduate students to grasp proof of the fact that 
the magnetic field outside the hollow superconducting sphere (superconducting shell) coincides with the field of a point 
magnetic dipole both when an uniform external magnetic field is applied as when a ferromagnetic sphere is located at 
the shell’s geometric center. The proof is based on the London’s equation transformation to an integral equation for the 
vector potential of system. The integral equation may be solved completely in case of the solid sphere. 
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Introduction 

The problem of uniform magnetic field expulsion, it’s value B ext is less then the critical one (B cr ), 
from the superconducting solid sphere volume (figure 1) was discussed in textbooks on the 
electrodynamics and on the solid-state physics many times. For a student - physicist it is the useful 
illustration of Meissner effect in elementary functions and without any idealization of situation, 
which is occurred, for example, in the problem of “infinite” slab in a parallel field. For the case of 
approximation X —» 0, where X is the field penetration depth in a superconductor we find it’s 
presentation in the initial university physics course already (Sivouchine, 1983). By this author of 
the textbook makes educated guess, that in an external field the sphere behaves like a point 
magnetic dipole with an u nkn own magnetic moment 91. Using of such educated guess seems to be 
quite acceptable in the general physics course however it is hardly right in textbooks for 
undergraduate students. Meanwhile the prompt message that the supercurrent modulus j(R) angular 
dependence is 


j(R) = j o (R)sm0, (1) 

(. R , 6 are the coordinates shown in figure l, jo is a function of only the R coordinate) may be 
occured, for example, in (Batygin and Toptygin, 1978) although this book is intended for a more 
advanced student. In this textbook the superconducting solid sphere is considered out of the 
approximation X -* 0 already, but in the situation when the London’s equation 

rot j = ~ f~72 B ( 2 ) 

AjtA, 
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Fig. 1 . Solid superconducting sphere of radius a in an external uniform magnetic field. The arrows are the 
total magnetic field B induction vectors in different points of the sphere’s equatorial plane 

is applicable (B(R) is the total magnetic induction vector, B = B cxt + B } , where B^(R) is the 

supercurrent field induction vector, c is the velocity of light). The method of solution in (Batygin 
and Toptygin, 1978) is the next: the projection of equation for supercurrent, 


rot rot j = -— j 




(3) 


on the e, basis vector of the figure 1 coordinate system, is, taking into account the supercurrent 
modulus independence on the azimuthal angle cp : 


1 . d 2 j 2 dj 1 d 2 j 1 cos 6 dj 1 . 

^ t 2 sin#^ dt 2 t dt t 1 dd 2 t 2 sin# dd t 2 sin# ^ 

(t = R/ A ). Equation (1) yields with (4) 

t 2 ' Jo + 2t • j' 0 - (t 2 + 2) • j 0 = 0. 


The general solution of (5) is (Nikiforov and Uvarov, 1988): 


/ 1 1 ) 

, , (1 1 \ 

1 

1 

J7o| 

e+k \E7) 


(4) 


(5) 


( 6 ) 


where k\ and kj are some constants being defined by boundary conditions. For the solid sphere we 
must demand j(R = 0) = 0. Expanding exponents in (6) into a series and inserting t = 0 there, we 
obtain k\ = ki = k. Further, using (2), student may calculate the total field induction value B(R) at R 
< a and, in particular, the value B(a - 0). Postulating that at R> a 


B,= 

B e = 



2-9f\ 



cos#, 


- B. 



■ sin 6 , 


(V) 

( 8 ) 
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and joining components of B at R = a - 0 and at R = a + 0 one can derive 

k = 3c- B ext a /16^>l 2 sinh( a/ A), 


( 9 ) 



( 10 ) 


and the problem is solved completely. 

The method of attack of Batygin and Toptygin was practically criticized in the article (Matute, 
1999). On author’s opinion, «the general formal solution is avoided)) in (Batygin and Toptygin, 
1978) and other textbooks, although «the problem can be solved by standard methods which is 
actually presented in the textbook)). In the article (Matute, 1999) the field equation 



( 11 ) 


is considered together with the Maxwell equations. Combining it and taking into account the 
symmetry of problem one come to some equation for the B r field component only. Solution of this 
equation is searched as a series in the Legendre’s polynomials P n (cosd). The B, component of 
vector B is searched as a series in terms of the functions P l n = dP n (cosd )/dd . 

The Maxwell equations solution for R > a is kn own (Jackson, 1962): 


71 = 00 

B, = B e« ■ cos e + J P n (cos d ), 


( 12 ) 



(13) 


where c n are still unknown coefficients. Equating (12) - (13) and the solution of (11) at r = a , they 
derive c n , by this c„ = 0 for all indexes n > 1. That is to say that when R > a, the field B j of 
superconducting solid sphere coincide with the field of a point dipole. Thus the fonnulated problem 
may be solved without any prior guess. The results of article (Matute, 1999) gave the groundwork, 
which is necessary for the exhaustive understanding of the solution of problem by an 
undergraduate student. However attempts to adapt the algorithm, having been developed in this 
work to analogous problem concerning the hollow sphere show noticeable increasing of the routine 
calculations volume, which are necessary for the proof of the super-current field B\ dipolar 
character. This proof became even more cumbersome when the magnetic dipole is added into the 
sphere’s cavity. Can an educator recommend to an undergraduate student the other way of solution, 
which will be not more cumbersome and will be interesting for the intellect? 

Hollow superconducting sphere (superconducting shell) in an uniform magnetic field. 

Let us consider the hollow superconducting sphere of inner radius h, outer radius a and volume V 
(figure 2). Let us prove the next prime statement: in the external uniform field B ext 
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Fig. 2. Hollow superconducting sphere of inner radius h and outer radius a in an uniform magnetic field. 
The context of arrows is the same as in the previous figure. 

the field of supercurrent BfR > a) coincides with the field of a point dipole. For this purpose we 
introduce vector potentials A i , A ext , A accordingly the fonnulae B } = rot/1, , B ext = rot A ext , B = 
rot 4. Owing to the cylindrical symmetry of problem we can calibrate it so that all of the three 
vector potential modulus are not depend on the azimuthal angle cp of our coordinate system. Then 
all of the three vector potentials divergences are zero in all space. According to (Batygin and 
Toptygin, 1978) in such a situation 

A(R) = A ext (R) + A j (R). (14) 

By this (Jackson, 1962; Batygin and Toptygin, 1978), 

w-Uffif- <15> 

As d\\A = 0, one can use the London’s equation in the fonn 

(16) 

4 Jtf 

Combining (14) - (16), we derive ( h<R< a): 


A{R) = A ext ( R) - X — ■ f^ldV , (17) 

ext 4^rzL 2 

where 

A„(R)-^f^sm0e r . (18) 

Let us confine ourself to R radius-vector, lying in the plane of figures 1-2. Projecting the equation 
(17) on the axis, which lies perpendicularly to this plane (Jackson, 1962), we derive the equation 
for the vector potential modulus already 
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A(R) = 


B mt • R . 


-sin#- 


A(R') 


6 -— • f ~r~ -: • cos cpdV'. 

AjtA 2 J. R-R 


(19) 


In dimensionless variables r = R / a, A = 


B e«a 


■ A , this equation may be rewritten as 



( 20 ) 


Here 3 = a 1 1 Ajt?i , V is the volume of unit outer radius hollow sphere in an imaginary space; it’s 
inner radius is h / a. We have the so named 2 - type nonhomogeneous Fredholm equation. One of 
theorems of the integral equations theory states (Korn and Korn, 1968): if the integral 

Fir) -v — r') dV’ is bounded function in V , then the number "E, exists such that at all 6 < § , 

v 

Neumann’s step-by-step approximations series 


A [l] (r) = r •sin#, 


j[n+ 1 ] 



•cos cpdV', 


( 21 ) 

( 22 ) 


is converging to a solution of (20) uniformly inside V. One can find in (Prudnikov, Brychkov and 
Marichev, 1986): 





12 


(23) 


and the boundedness of both function G{r ) as a fortiori function F(r) < G(r) at 0 < r < 1 became 
evident. Suppose the formula A = A(0)/ -\\ -(T!T cr ) 4 is correct (Ashkroft and Mermin, 1976) where 

r cr is the superconducting transition temperature. We conclude that some temperature interval T x < 
T < T cr exists where the series of functions (21) - (22) is converging to a solution of equation (20). 
Near the T cr , where 6 -> 0: 


r ' sin# 


A(r) ~ A [2] (r ) = A m (r) + A A(r) = r ■ sin 9 - 3 • f—- — • cos cpdV' 

d \r-r\ 


(24) 


Using the expansion of l/|r - r I in spherical harmonics (Jackson, 1962) we obtain 


oo / 


1 


AA(r) = -Ajz3 • y ^j~[ Yl ^ 6,(p = °) ‘ / F ii r ’ r 'y adr '$ Y Lmi 6 ' ’V}' sin & cos cp'dS' (25) 

/=0 m = -I r'=h/a S 


, where S is the unit sphere surface, F { (r,r') = min(r,r')' /max(r,r') + . As 
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Lm (0,(p = 0) ■ fY l m (8' ,(p')smd' ■ cosqp'dS' = 6 lx -sind, 

m=—l s 

what is verified easy (d m> n is the Kronecker symbol), then 

A [2] (r) = f 2 (r)-sin 8, 


(26) 


(27) 


where f 2 (r) = r ——d • J'F l (r,r')r 13 dr' is a function only of the r variable (there is no necessity 

^ r'=h! a 

to calculate / 2 (r) in an explicit form although it may be made very easy). It follows from (22), (26) 
and (27) that formula with the structure analogical to (27) will be reproduced at any step of iteration 
process with character replacement 2 -> n in (27). So, when the superconductor is cooled down 
from T cr temperature, only the radial distribution of field is changed but it’s angular distribution is 
unchanged. We conclude that in all the region of Neumann’s series (20) - (21) convergence, at h < 
R <a, 


A(R) = f(R)-sinS-e r 


(28) 


where J{R) is a function only of the R variable. Let us insert (28) and (16) in (15) at R > a. It gives: 


A-(R > a) = - 


B ca a 

R 2 • 6A 2 


\fmr 


dl I • sin 8 ■ e 


(29) 


This is the vector potential of a point magnetic dipole. We have proved the prime statement at the 
temperatures T x < T < T cr . That will do as we now have a right to use the function (6) in this 
temperature interval and further to confirm that formula (6) is correct at an arbitrary temperature. 
Now if a student knows indefinite integrals of the elementary functions, he will calculate the value 
A(R) at h < R < a explicitly without problems. For this purpose it’s necessary to insert (1), (6) and 
(16) in the both parts of (20) and to equate coefficients at the same degrees of r. We present the 
final result: 


f(R) = 

25 


(k A 2 ) „ ^ 


(k k 2 \ m 


R R 2 


U-e 


- 1 - 7 

R R 2 


W-e 


h 1 h h 2 h 

where U = — T + — + 1, W = — T -+ 1 and 

3k k 3k 2 k 


a-h h-a 

S = U-e T -W-e^ 


(30) 


(31) 


Inside the cavity, the formulae (15) - (16) give: 


R R<h) -^ 


! a 


ff(3)dt • 


sxnd-e„ 


\ h 


(32) 
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and after a simple calculations 


B(R<h) = 2a-B ext /AS. 


( 33 ) 


Expanding exponents in (31) as a series in the parameter e = a / X, one can prove that at any h and a 
it will be B(R < h) < B ext always. Hollow superconducting sphere in a field of point dipole located 
at it’s center and the generalization of problems. 


In the article (Hurault and Pincus, 1969), the problem was formulated inter alia: to detennine 
magnetic field inside the infinite superconductor surrounding a monodomain ferromagnetic sphere 
of radius / and of magnetization M s . The solution of equations (3) and (16) for this 



Fig. 3. Hollow superconducting sphere of inner radius h and outer radius a in the field of a magnetized 
sphere with the dipolar moment p m . The context of arrows is the same as in the previous figures. 

problem, resulting in this article is 


A(R > I) = 4 jM s 


l 3 R/A + l 

R 2 3 + 3UA + I 2 /A 2 


e -(R-t)U 


fmd-e v . 


(34) 


The method being developed in the previous paragraph allows to provide the necessary background 
for this fonnula which is absent in (Hurault and Pincus, 1969). We will make it for the more 
common situation of the superconducting shell in the field of a ferromagnetic sphere with the 
magnetic moment p m located in it’s center (figure 3). Just the same method as in the previous 
paragraph leads us to the equation, 


A(R) = sin 6 ■ e v - X — ■ C^£±dV' 

R 2 v 4 jtA 2 J. R-R 


(35) 


instead of the equation (19). The dimensionless variables are here r = R/a, A = a 2 ‘Alp m . Then 


A(r) = -^-sin 6-d - r ^( r ) - cos cpdV'. 

r Jl r — r 


(36) 
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Further consideration duplicates the text beginning from formula (22) to formula (28) with the 
rewritting of function A ll] (r ) = r“ 2 -sin 6 instead of function A m (r) = r -sin 6* overall. Inserting (1), 
(6) and (16) to the both parts of (36) and equating coefficients of the same degrees of r, we obtain 


f(R) = 


Pm 

a 2 -s 


7 A 

Ac \ 

R-a 

(A Ac \ 

a-R - 


-* 1 ) 

■e A + 

U + * 2 J 

•e A 


When a —* oo and £ = h we reproduce the result of (Hurault and Pincus, 1969). 
The field inside the cavity (/ < R < h): 


Pm * sinhl 


B(R) = B lm (*)-- 


' a-h' 


hA 2 S 


(37) 


(38) 


where Bf m (R ) is the field of ferromagnet outside it. 

Let^4i(/f) is the solution of the equation (17) for the superconducting shell in the field of the 
uniform magnetic field B ext and A 2 (R) is the one’s solution for just the same shell in the field of 
point dipol being located in the center of shell (B\ = rot4i , Bi = rot42 , B = B\ + BA). For the case 
of very weak field B (B < B cr is elsewhere), A(R) = A\(R) + A2(R) will be the vector-potential for 
just the same superconductor in the field both of that dipol inside it and of B cxt outside it. A 
consequence of the results stated above is the fact that the formula (28) will be correct also in this 
most common situation. The detailed calculations of the field distribution may be easy executed 
with the help of formulae (30) and (37). 

Conclusions 


Using results of the integral equations theory allows to put the novel method of solution of the 
problem about the superconducting sphere in a magnetic field into pedagogical practice. This 
method is alternative to the traditional one, when the London’s equation solution inside a 
superconductor and the MaxwelFs equations solution outside it are joined on it’s surface. The 
advantage of this method over the tradition one lies in the fact that the change from the problem 
about the solid sphere in an uniform field (problem N°l) to the problem about the hollow sphere in 
an uniform field plus the field of a dipole (problem N°2) doesn’t bring the solution’s significant 
amplification. One may say that now not only the problem N° 1 is open to understanding of a student 
but also the problem N°2. Nevertheless if an educator on the lecture wants to confine oneself to the 
situation of the solid sphere, he can now demonstrate the solution of equation (20) which is given 
in Appendix instead of to integrate the contents of (Matute, 1999) and (Batygin and Toptygin, 
1978). 
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Appendix 


When h = 0, the integral equation (20) is solved very simply and there is no necessity to cite the 
fonnula (6). For this purpose let us insert (28) in (19). In combination with (26) it gives the next 
integral equation for the function a)(r) = 2/ (r)! B ext a: 


co(r ) = r - — • J — • J aft) • fdt + r J co(t)dt 


(39) 


As aforesaid, the functional sequence eo [n \r), where A [n] (r) = co [n] {r) ■ sin# convergents to a 

solution of (39) at 0 < r < 1 and at an infinitesimal e. Being made some iterations, student may 
notice that 


c o [ "\r ) = 


i=n 



f^A-' 


• e 


2(i-l) 


(40) 


where b n jj are some coefficients. Formula (40) is proved elementary by the method of mathematical 
induction. Suppose off) is infinitely differentiable function at r = 0 both on r and on e variables. So 
when n —> oo, coefficients of the series (40) must tend to the coefficients 
(1/ k\-m\) ■ d k+m to / d k e d'"r of the function afr) Taylor’s series. It signifies that we may 

£=r =0 

change the sequential order of indexes in (40): 

j=n { i=n 

“'"’M - %l b P 

j =1 \ i=j 


2(/-l) 


• r 


2y-l 


(41) 


So we may search the solution of (39) as 


J-™ 

(o{r) = ^bj- r 2hX 


(42) 


i=n 

( bj = lim,,^ ^b" i £ 2(l ^ ]) )• Inserting (42) in both parts of (39) and equating coefficients at the same 

i=j 


degrees we derive 
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p 2 J=cc h 

6,-1 + —-t —= 0, 

3 falj 


b i+ 1=--- -b r 

J 2y(27 + 3) ' 


Now let’s insert (44) in (43). It gives 

2 


b , 


1 H-1-1-h ... 

2-3 2 ■ 4 • 3•5 2 - 4 - 6 - 3 - 5-7 


= 1 , 


(43) 

(44) 

(45) 


so that b]_ = 81 sinlu: and consequently 


( 2 3 4 5 6 7 

£ T £ T £ T 

r +-+-+-+... 

2-5 2 • 4 ■ 5■7 2 • 4 • 6 • 5 • 7•9 , 

37?, / 1 4 3 6 . . 5 8 . 7 \ 3 / cosh(£ • r) sinh(£ • r) \ 

= — L - -E-r + —{e-r) +—(e-r) +—{e-r) +... = •' 


£ (3 5! 

Definitely 


7! 


9! 


sinhf 


£ • r 


( £ ' r ) 


• (46) 


/ 


^4(7? < a) = 


3^ ex t« 


2 sinh(a / A) 


cosh(7?/A) sinh(/?/ A) 


7?/A 


(7?/A) 2 


sin6»-^. 


(47) 


Inserting (46) to (29) we immediately come to (10). 
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